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i-£h Abstract: We consider a set of scalar fields, consisting of a single flat direction and one or 

several non-flat directions. We take our cue from the MSSM, considering separately D-flat 
2 and F-flat directions, but our results apply to any supersymmetric scenario containing flat 

directions. We study the field fluctuations during pure de Sitter inflation, following the 

i_i evolution of the infrared modes by numerically solving the appropriate Langevin equations. 

| We demonstrate that for the Standard Model U(l)y, SU(2)l or SU(3) C gauge couplings, 

as well as for large enough Yukawa couplings, the fluctuations along the non-flat directions 

^ effectively block the fluctuations along the flat directions. The usual expected behaviour 

(4> 2 ) oc N, with N the number of efolds, may be strongly violated, depending on the 
coupling strengths. As a consequence, those cosmological considerations, which are derived 
assuming that during inflation flat directions fluctuate freely, should be revised. 
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1. Introduction 

The scalar potential of a supersymmetric gauge theory has always a number of flat direc- 
tions, along which the potential vanishes identically. In particular, the scalar field space 
of the MSSM (Minimal Supersymmetric Standard Model) is 49-dimensional, and there are 
some three hundred flat directions; these have all been classified (for a review, see [1]). The 
flat directions are rays in field space, restricted to lie in some subspaces. More concretely, 
flat directions correspond to configurations where some of the field values are related to 
each other while the rest are set to zero. One consequence is that most of the flat direc- 
tions are not simultaneously flat; if there happens a field fluctuation along some given flat 
direction, many other potentially flat directions will no longer remain flat. This is true in 
general for any supersymmetric gauge theory with scalars. 

In the absence of supersymmetry breaking, flat directions are perturbatively safe as 
they are protected by a non-renormalization theorem [2]. However, they can be (and are) 
lifted by supersymmetric non-renormalizable terms. For each flat direction, the lowest order 
non-renormalizable operators lifting the flatness are known and classified. Supersymmetry 
breaking will also induce non-flatness through the mass terms and A-terms, but for the 
purpose of the present paper they are not important. 

For concreteness, we will focus on the MSSM, although everything that we have to say 
will generically hold for any supersymmetric model. A central ingredient in the cosmological 
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considerations of flat directions is often the tacit assumption that during inflation, fields 
can fluctuate along a given flat direction essentially like massless fields in de Sitter space. 
The variance of a massless field (p m pure de Sitter space is given by [3] 

<^ 2 > = 2r* H ° N > (L1) 

where N is the number of e-folds and H Q is the Hubble parameter during inflation. Fluc- 
tuations are assumed to continue growing in amplitude until the non-renormalizable terms 
take over, making the effective mass V" > H 2 , and hence stopping the mean-square- 
root-amplitude from evolving any further. Quantitatively, for a flat direction with mass 
m <C H , the asymptotic behaviour t — > oo of the variance is given by [4] 

^ 'l-^f-^)) , (1-2) 



8ir 2 m 2 V \ 3H Q 

Hence if m ~ H Q , the mean-root-square amplitude of the fluctuations would saturate to 
a constant value. Assuming that the effective mass is only due to the non-renormalizable 
term(s), one would obtain a very large flat direction field amplitude that depends on the 
dimension of the non-renormalizable term(s). Taking V ~ (p 6 /Mp as an example, with 
M p ~ 2.3 • 10 18 GeV the reduced Planck mass, once the condition V" ~ H 2 is attained, 
one would find a final amplitude of (poo ~ (H Mp) 1 / 2 3> H Q . Note that this additionally 
assumes that inflation lasts long enough as reaching the asymptotia in this particular case 
would require N ~ M p /H Q efolds. However, this is not quite true since from (1.2) we 
see that equilibrium is attained when V ~ H%, corresponding to cp CCL ~ (H^M 2 ) 1 / 6 = 
(Ho/Mp) 1 / 6 cftoo. The equilibrium probability distribution V is not peaked at <p cq but is flat 
with V ~ exp(-F/i?£) [5]. 

However, during inflation all fields fluctuate, including the non-flat directions. Since 
the existence of a given flat direction is conditional upon most of the fields staying at the 
origin, one may wonder how likely it is that a large amplitude along a flat direction can be 
obtained, considering that the flat direction is actually coupled to non-flat directions, which 
are also fluctuating. They may in effect provide an effective mass for the flat directions, 
thus preventing the spreading of the flat direction variances. We will demonstrate that in 
most cases, depending on the gauge and Yukawa coupling strengths, the small fluctuations 
of the non-flat directions are indeed sufficient to block the growth of the mean-root-square 
amplitude along the flat direction. 

In the present paper we address this issue by separating the ultraviolet (UV) and 
infrared (IR) field modes and solving the appropriate infrared Langevin equations for a 
restricted set of MSSM fields in de Sitter space. We consider the components of the Higgs 
fields H u and a slepton field L that are mostly non-flat field directions but, in a certain 
combination, also include a flat direction; we also include other non-flat field directions. 
Our conclusions will be, in any case, extrapolable to other scenarios with flat and non-flat 
directions interacting. 

The contents of the paper are as follows. After the introduction presented here in 
Section 1, the fields and the potential are written out in Section 2. In Section 3 we review 
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the formalism of the Langevin equations and the separation of the IR and UV modes, 
and describe our procedure for the numerical implementation. In Section 4 we present the 
outcome of the numerical solutions we obtain, and discuss the role of the coupling strengths. 
We treat F-flatness and D-flatness separately. In Section 5 we draw our conclusions. 

2. Coupling of flat and non-flat directions 
2.1 The scalar potential 

A supersymmetric potential is given by the sum of the F-term and the D-term: 

2 

4- 

2 



2 

Vs = Vp + Vd = £ ® + l Y^9 2 a D a D a , (2.1) 



a 



where (pi stands for all the scalar fields. In the MSSM the superpotential reads 

/ = \ u QH u u + X d QH d d + X e LH d e + fiH u H d , (2.2) 

where Q and L are respectively the squark and slepton doublets, and H u and H d are the 
Higgses. Here {g a } are the Standard Model gauge couplings, {Aj} the Yukawas, and 



D a = J24 Ta ^i (2-3) 



is the D-term. If for some set of field values F = (D = 0), we call the corresponding 
direction in field space F-flat (D-flat). 

During inflation, the fields will typically obtain an induced mass term with m ~ H Q . 
However, the magnitude and also the existence of such a term depends on the details of 
the inflationary sector (for a review, see e.g. [1]). For instance, in D-term inflation the 
Hubble-induced mass term is absent. In what follows we will ignore it, which in the light of 
our results is the conservative assumption since the Hubble-induced mass will usually only 
add to the blocking of fluctuations (although it is also possible that the Hubble-induced 
mass term has a negative sign). 

To be concrete, we will focus on the leptonic sector of MSSM only and take our cue 
from the simple LH U flat direction, given by 



h u = : ,l= : (2.4) 




with all the other scalar fields = 0. It is easy to verify that the configuration (2.4) is both 
F-flat and D-flat. Thus we will assume that the squark fields do not fluctuate; since they 
are not coupled to the leptonic sector, their fluctuations would affect only the Higgses (but 
nevertheless also the LH U flat direction). In addition, we also set H£ = for simplicity. 
After an SU{2) rotation, the remaining degrees of freedom are: 



H d = n > L =\ h e , 0, (2.5) 
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where we have denoted the prospective flat direction as = (j). As we will see, some of 
these dofs are, qualitatively speaking, dynamically redundant. 
Given the fields (2.5), we find 

V F = \ 2 e (\eH°\ 2 + \ee\ 2 + \eH°\ 2 + \ve\ 2 ) , (2.6) 

while 

^ 2 = ^(M 2 -M 2 + |tfT-M 2 ) , 

Di = \ (-M 2 - |e| 2 - |#°| 2 + 2|ef + \^\ 2 ) (2.7) 

are respectively the SU{2)l an d U(V)y D-terms. As is conventional, we neglect the //-term 
which is of the order of susy breaking mass. 

From (2.5) and (2.7) we readily observe that the flat direction is given by v = (/), 
provided the right-handed slepton e = 0. 

2.2 Fluctuations breaking F-flatness 

Note that in general the flat direction (f> mixes with the non-flat part of v in the kinetic term 
(a similar mixing arises also through susy breaking mass terms). This is just an artifact 
of the chosen notation which presents a computational nuisance, since in the evolution 
equations for fluctuations, one would have to diagonalize the kinetic term at each time 
step. However, this can be avoided if we choose <fi to be real (with cj) — > (fr/Vty and write 

„ = -L(0 + ifr) , (2.8) 

where h is real and represents the difference v — H®. Then the kinetic terms of all fields will 
remain diagonal. We do not believe that the restriction (2.8), which represents a subclass 
of all possible fluctuations, changes the outcome of the dynamics of the F-term in any 
qualitative way. However, in the D-term the choice (2.8) decouples the flat direction </> 
from the non-flat directions completely. Therefore the effect of the fluctuations of non-flat 
directions on the flat direction in the D-term must be considered separately and will be 
addressed below. 

Let us now inspect (2.6) and (2.7) to see what sort of generic couplings we obtain. 
Let us also denote the non-flat directions generically by Xi- We see that from the F-term 
one obtains potential terms of the type |xiXj| 2 with i 7^ j, as well as a coupling of the flat 
direction to non-flat direction e through \ve\ 2 = (eft 2 + h 2 )\e\ 2 . With the form (2.8), there is 
no ^-dependence in the D-terms, and one obtains potential terms of the form \xiXj\ 2 with 
i 7^ j, as well as terms like \xi\ ■ Therefore, the generic form of the potential when some 
non-flat excitations are coupled to the flat direction in the F-term, can be summarized by 
the example 

Vs = \\ 2 M 2 + h 2 )e 2 + \&h* + \ 9 l^ + 4e 4 - 4h 2 e 2 ) + ^ , (2.9) 
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where we have chosen the fields to be real (and properly normalized) and have set e = 
since their contribution to the potential is qualitatively similar to e; as such, they 
just represent some additional fluctuating dots that couple to each other but not to the flat 
direction. The last term in (2.9) is the non-renormalizable term appropriate to the LH U flat 
direction, which is known to be lifted by d = 4 operator in the superpotential. For other 
flat directions the dimension can be different but always d > 4; all flat directions are lifted 
by operators with d < 9 [6, 7]. It will turn out that, for most cases, the non-renormalizable 
term is irrelevant for determining the final amplitude along the flat direction. 

We will treat a potential of the type (2.9) as a generic example describing a flat 
direction <f> and two non-flat directions that may all fluctuate during inflation; however, 
as noted above, the full spectrum of fluctuations is not included as the flat and non-flat 
directions remain decoupled in the D-term. In what follows, we assume that the Hubble 
rate during inflation is much larger than the susy breaking scale. Hence we will neglect 
both the supersymmetry breaking mass terms and A-terms in the potential. 

2.3 Fluctuations breaking D-flatness 

To study the effect of fluctuations inside the D-term, it is not convenient to consider the 
flat direction as a background solution. Rather, we simply decouple the F- and D-terms 
by setting the Yukawas to zero and consider the field fluctuations in the potential 

V D = \g\L*-H*)\ (2.10) 

where L and H have been assumed to be real dofs for simplicity. Their kinetic terms are 
diagonal. The flat direction is the ray L = H, and we will consider the evolution of the 
probability distribution due to inflationary fluctuations in the (H , L) plane. The fluctua- 
tions of other (non-flat) fields are ignored for simplicity. Although (2.10) is motivated by 
MSSM, we will treat the coupling g as essentially a free parameter, focusing however on 
values that are appropriate for the MSSM. 



3. Stochastic dynamics of fields during inflation. 

Let us now study the field fluctuations during inflation, assuming a constant Hubble rate 
H Q . Any field in de Sitter (or quasi-de Sitter) can be decomposed into IR modes (with 
momenta k < eaH Q ) and UV modes (with k > eaH D ), where e is a constant smaller 
than unity, e < 1. The IR modes have non-trivial stochastic dynamics [8] due to the 
continuous influx of UV modes around the Hubble radius, which can be considered as 
classical stochastic fluctuations as long as e is sufficiently small; any value e < 1 will do. 
In this manner the UV modes impart continuous random 'kicks' on the IR sector that 
are superimposed over the usual deterministic evolution. The dynamics are described by 
multi-field Langevin type equations, from which one can extract the probability distribution 
function (PDF) of the fields' fluctuations. 
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3.1 The Langevin equations 

We begin by briefly reviewing the formalism of the IR stochastic dynamics during infla- 
tion [8] (see also eg [9, 10] for some early references on the subject). Let us define the IR 
and UV parts of a scalar field <p(x, t) as 



d 3 k 

(2tt) j 
d 3 k 
(2tt) j 



IR (x,t) = / ^e-^cMWM (3.1) 
JV (x,t) = / -_^e- kx ^ k (t)[l-W(k,0], (3.2) 



where the function W^(k, rj) is an IR-filter (window function) that is subject to constraints 

<-l, (3.3) 



w(Kt) |k|<<Q(f) > i, w(Kt) |k|>>Q( V 

dk 



with Q(t) a time-dependent scale. For inflation, the relevant characteristic scale is the 
(comoving) Hubble radius, H = aH Q , which provides a natural border between the UV 
and IR dynamics of the modes of a light scalar field. We will thus identify Q(t) = eH, with 
e < 1 . Since the energy density of the fields will be sub-dominant they will not back-react 
on the metric during inflation, and we will ignore gravitational perturbations; for their 
inclusion in a stochastic formalism see [11, 12]. 

Consider now a set of interacting fields in de Sitter space, {fa}. Each field can be 
decomposed as <^>i(x, t) = $j(x, rj) + t), with $j and fa the IR and UV parts obtained 
according to eqs. (3.1)-(3.2). In order to most accurately follow the dynamics of the IR 
dof, $j and should be considered as independent variables (though linked through the 
eom) and a Hamiltonian formulation is the most natural framework to use. Calling 7Tj the 
conjugated momentum of fa, the eom in the hamiltonian picture are 

fa = 7Ti, 7Ti + 3H -Ki = \\7 2 fa - diV , (3.4) 

cr 

with diY = dV/dfa. We then IR/UV decompose fa and Hi independently as fa(x,rj) = 
$j(x, rj) + ¥?«(x, rf) and 7Tj(x, rj) = IIj(x, rj) + <5-7Tj(x, rj), and introduce such a decomposition 
into eqs. (3.4). We provide more details in the Appendix. One finally finds that the 
dynamical equations for the IR dofs are 

= n i + s 4 W (x,? ? ), (3.5) 
Hi = -3H Ui - DiV + s\ n) (x, r?) , (3.6) 

with Di = d/d^i, V the potential taken as a function only of the IR components, i.e. 
V = V({$j}), and 

#W) = J ^- 3 e-^fa(k,t)W(k,t) (3.7) 



* l (M) = ' j^y 3 e-^n t .(Kt)W(k,t). (3.8) 
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Note that since we are interested in the leading order IR behavior we have dropped the 
gradient terms from (3.5) and (3.6), and have also ignored any corrections coming from 
the IR/UV decomposition of an interacting potential V. 

Eqs. (3.5) and (3.6) are operator equations. The terms e _lkx ^fc, e _ * kx 7Tfc should be 
understood as afce~ lkx /^(k) + a|.e +lkx /^(k), a/ c e _lkx / 7r (k) + a^,e +jkx /*(k), respectively, 
with afc, a\ the usual creation/annihilation operators, and /</>(k), /^(k) the field mode func- 
tions. However, if the window function is chosen appropriately such that the UV/IR split 
occurs shortly after the relevant modes have crossed the horizon, the phase of the mode 
functions entering (3.7) and (3.8) becomes almost constant and then and com- 
mute with each other at different times. They can be considered as classical stochastic 
forces and the equations for the IR dofs are therefore Langevin-type equations. Ignoring 
the self-interactions of the UV modes makes the stochastic terms gaussian random fields 
so that all the statistical information is encoded in the correlators 



;;>(t^) = (a<«(x,t)*f (x,0) , (3-9) 



ij 

IJ 



;,V0 = (^W)4"W)) • (3-10) 



Note that since we are not interested in the spatial correlations, we are only considering the 
correlators at the same spatial point, which in reality corresponds to a region of physical 
volume V ~ V#o- P° m ts separated by physical distances L > l/H Q are essentially 
uncorrelated. 

Apart from the classical evolution dictated by the deterministic parts of (3.5) and 
(3.6), the stochastic forces acting over a small time interval 5t <C 1/H Q displace the fields 
by 

rt+St 

S] 



^ ] (t')dt', (3.11) 



t+st 

(tt) 



Sin = / sf\t')dt'. (3.12) 



Therefore, the correlators we really need are 

rt+St rf+5t 



S^(t,t';St) = J ^ sf(r,r')dTdT', (3.13) 
f t+8t r t'+St 

S*\t,t>;St)= / s^(r,r')drdr'. (3.14) 



J t JV 

In pure inflationary de Sitter background, the solution to the mode equations (with bound- 
ary conditions matching Minkowski modes at k — > oo) are well known, see for instance [13]. 
Using such mode functions, /^(k, t) and f Wi (k,t), evaluated at k = eaH Q , and choosing a 
step-function for the IR-filter, W(k,t) = 0{eaH o — k), the equal-time correlators we need 
are found 1 to be 

Sif(t, dt) = Sa (1 + 6 3 )|| dt , S\f{t, dt) = Sa ^^ydt, (3.15) 



1 See the Appendix for more details. 
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valid only in the limit dt <C 1/H Q . The choice e < 1 then ensures that the classical inter- 
pretation of the stochastic noise terms is valid. We also see that for e < 1 the momentum 
noise term is suppressed and the main stochastic component lies in the fluctuations of 4>. 

3.2 Numerical implementation 

Rescaling the field variables as 

->■ $i = ^/Ho , S W = S W/ij2 

n» -> fii = Ik/ Hi , sW -> sW = sW/F„ 3 

^({^})^^({^})/^o 4 , (3-17) 
makes it possible to express the Langevin Eqs. as depending only on dimensionless variables, 

^(N) = U i (N) + sf, (3.18) 
fL'i(N) = -SHofliiN) - D#({%(iV)}) + s\ w) , (3-19) 

with ' standing for derivatives with respect the number of e- folds N = J Hdt. 

In order to solve these equations in a computer, we just need to discretize them by 
choosing a small time step, i.e. dN <C 1. Dropping the tildes for clarity of the notation, 
one arrives at the discretized (dimensionless) iterative equations 

$i(N + dN) = $i(N) + Ui(N)dN + (3.20) 

Ui(N + dN) = Ui(N) - 3Ui(N)dN - A^({^(A r )})rfA r + sf , (3.21) 

where the stochastic terms are drawn from a gaussian random distribution with correlators 



sif(dN) = (S^sf) = 5 tJ (1 + e 3 )^ , (3.22) 



s\]\dN)^{S^sf) = 5 l3 e^ 2 . (3.23) 

Eqs. (3.20)-(3.23) characterize completely the dynamics of the IR dof of any set of 

scalar fields during de Sitter inflation. Note that the only scale of the problem 2 , the 

inflationary Hubble constant H Q , has been scaled out, so Eqs. (3.20)-(3.23) are indeed 

scale-free equations. Any dimension-full functional built by powers of fields, for instance 

(f) 2 , will then be measuring an amplitude in units corresponding to the same powers of 

H Q . Eqs. (3.20)-(3.23) are therefore universal in this sense, since the physics they describe 

is independent of the inflationary scale. It is also worth mentioning that the stochastic 

(x) 

character is time-independent, since the correlators S 1 ^- do not depend on the 'time' N, 
but only on the step dN. The dynamical behaviour of the IR dof should, of course, not 
depend on such step dN. Thus, after solving the system for a given step dN <S 1, one 
should always make sure that the same dynamics is recovered by choosing, for instance, 
dN/10. If that is not the case, one must then decrease further the step dN, until finding 
that the fields' statistical properties are insensitive to further decrements. Checks of this 
nature have been performed on all the numerical results we present in the next section. 

2 Of course, non-renormalizable terms in the potential will also introduce new scales, but as said before, 
and as we demonstrate with the numerics, such scales never play a role in the problem under study. 
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Figure 1: Time evolution of the Probability distribution function V(<f) of the Flat direction 
<j>. The PDF at each time/e-fold N was reconstructed with 10 5 independent runs, all for the 
case A = gi = g<i = 1. Each PDF shown in the sequence was obtained by dividing the unit- 
normalized histogrammatic distributions of <j> by its field bin width at each time. The colors 
represent different stages of the evolution. Purple-to-DarkBlue, from N — 0.1 to N = 1 e-folds; 
DarkBlue-to-LightBlue, from N = 1 to N = 3.5 e-folds; Green, from N = 3.5 to N = 10 e-folds; 
Green-to-yellow, from N = 10 to N = 100 e-folds; and Yellow-to-Red, from N = 100 to N = 1500 
c-folds. 



4. Numerical results for fluctuations 
4.1 Breaking F-flatness 

We have solved the previously discretized equations for the potential (2.9). In what follows, 
we will denote the flat direction as (ft while Xi?X2 are the non-flat directions. We have chosen 
for the initial conditions 4> = xi = X2 = at N = 0. For a given choice of the couplings 
{A, g±, 52}) we have made n r (^> 1) independent runs with the same initial conditions. Thus 
each run represents a different realization of the time evolution of the fields' fluctuations. 
Collecting at a given moment the amplitude of the fields from all the independent runs, 
we can then reconstruct the probability distribution function (PDF) V{<p) of any field 99 at 
such moment. By definition Vi^p'^d^p represents the probability of the field (p to have an 
amplitude within [ip' , <p>' + dip]. Thus V(<p) can be obtained from the numerics by dividing 
the (unit-normalized) histogrammatic distribution of ip by its field bin width. Doing this 
at successive moments, we can track the time evolution of the PDF of each field. 

If there was no coupling between the flat and the non-flat directions, 4> would freely 
fluctuate with an increasing variance (4> 2 ) ~ NHl, until the non-renormalizable term 
</> 6 /Mp would become important and equilibrium would be reached. The non-flat directions 
also fluctuate, obtaining variances of the order ~ Q.IH 2 in a matter of few efolds (bottom 
lines in Figure 2). Because of the coupling of cfi to the fluctuations of non-flat directions 
Xi, the flat direction can then obtain an effective mass of order ~ 0.1XH o . Unless the 
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Figure 2: Evolution in time of the variance of the flat direction (ft in the strong coupling regime 
A = <7i = <72 = 1. Also shown is the evolution of the variance of the non-flat directions. The different 
plots of the flat direction correspond to different reconstructions of the evolution of (eft 2 ) versus N, 
as inferred from realizations with different number of independent runs, n r = 10 3 , 10 4 , 10 5 and fO 6 . 
Even for the poorly sampled case of n r = 10 3 , the physics is already well captured: the growth of 
the flat-direction variance stops after ~ 1000 efolds, reaching an asymptotic value of ((ft 2 ) s» 3.5H 2 , 
much smaller than the expected value NH 2 / Air 2 for a free flat direction massless field. 

couplings are very small, (ft can no longer be considered as an effectively massless field in 
de Sitter space. As a consequence, its fluctuations will be blocked. This is demonstrated 
in Figures 1 and 2, which depicts the case A = g\ = g% = 1. 

In Figure 1 we show the evolution in time of the PDF of the flat directions. In Figure 2 
we show the time evolution of the variance of each field. As can be clearly seen in the plots, 
it takes quite a long time, of the order of ~ 1000 e-folds, to reach a stationary regime. Once 
in this regime, the form of the flat direction PDF does not evolve in time anymore (red 
curves in Figure 1), and the variance growth stops, reaching a final asymptotic amplitude 
(see Figure 2). Note nevertheless, that during the first few tens of e-folds, the variance grows 
linearly with N, but still much more slowly than in the purely massless case. After the 
first hundred e-folds, the growth of the variance slows down and approaches the stationary 
regime. After N ~ 1000 efolds, the flat-direction finally enters into the stationary regime 
and its variance reaches an asymptotic constant amplitude ((ft 2 ) ~ 3.5H%, much smaller 
than the expected amplitude for a free massless field ((ft 2 ) = {N / Att 2 )H 2 . 

In all of this, the non-renormalizable term plays no role, since the flat direction field 
never manages to fluctuate further than few units of H Q . For A = 1, the amplitude of the 
flat direction in the stationary regime, is ~ 0(y/M p /H ) orders of magnitude smaller than 
the value needed for the non-renormalizable term to become relevant. 

However, all this depends on the magnitude of the couplings. In Figure 3, we show 
the time evolution of the variance of the flat direction variance with different values of the 
Yukawa coupling A, whilst g\ = <?2 = 1 as before. For A < O(10 -2 ), the flat-direction 
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Figure 3: Time evolution of Variance of the Flat direction, for different Yukawa couplings A, from 
A = 1CF 6 to A = 1. For the smallest cases, A < 0(1O~ 2 ), the flat-direction fluctuates as in the case 
of a massless free field in de Sitter space. For A ~ 0(1O -1 ), the flat-direction continues to grow for 
at least thousand e-folds, albeit more slowly than in the massless free case. For A ~ 0(0. 1) — 0(1), 
the freedom of the flat-direction to fluctuate is gradually blocked as A increases. 

appears to fluctuate as one would expect from a massless field in de Sitter space. For 
A ~ O(10 _1 ), however, we already see some deviation from the purely massless case, since 
the flat-direction continues to grow for at least few thousand e-folds, albeit more slowly. For 
values of A between 0.1 and 1, we observe gradual blocking of the flat-direction fluctuations. 

Thus, at the qualitative level, we may conclude that fluctuations within the F-term 
block the fluctuations along the flat directions that involve the third generation squarks and 
sleptons, while those involving only the first generation are likely to remain free to fluctuate. 
However, this is a conclusion based on the F-term alone and must be supplemented by the 
additional blocking provided by the D-term, as will be discussed in the next subsection. 

We should like to stress that the PDF properties do not depend on the richness of the 
statistics. As can be seen in Figure 2, n T = 10 3 runs capture already all the essentials. 
There we show the evolution of the variance also for n r = 10 4 , 10 5 or 10 6 independent 
runs, and all plots show the same behaviour, the only difference being that the richer 
the statistical sampling (i.e. the bigger n r ) the less noisy the plots are. As yet another 
check (not shown in the plots), we have compared the outputs from the same number 
of independent runs but different step size dN. We have found that if the simulation 
incorporates more than 100 stochastic kicks per e-fold (i.e. if dN < 0.01), the properties 
of the PDF remain the same as compared to the dN = 0.01 case, no matter how small we 
make dN. For simulations with n r = 10 5 or more independent runs, we cannot even detect 
any difference between choosing dN = 0.01 or dN = 0.001. 

4.2 Breaking D-flatness 

Let us now focus on fluctuations within the D-terms. We set the F-term (i.e. the Yukawas) 
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Figure 4: Time evolution of the variance of the field L {H is practically identical), for different 
gauge couplings g in the potential (2.10). The larger the coupling is. the stronger the blocking of 
the flat direction fluctuations, and thus, the smaller the amplitude reached. With relatively strong 
SU (3) interactions, the fields reach a stationary regime after only N ~ 200 — 300 cfolds. In the 
SU (2) case, the fields are just entering into the stationary phase after the first ~ 1000 — 2000 cfolds, 
whereas in the U(l) case, the amplitude of the fields is still growing after N = 2000 cfolds, albeit 
with a much smaller amplitude than in the free massless case (straight line in the Figure). 

to zero and adopt the potential (2.10). We have thus two fields, H and L, interacting with 
the coupling g, for which the configurations H = ±L represent a flat direction. 

In Figure 4 we show the evolution of the variance of L and H for different strengths of 
the coupling between them. As in the case with an F-flat direction, we see again a gradual 
blocking in time of the freedom of these fields to fluctuate. The fluctuations of H and 
L are restricted more severely the stronger the gauge coupling. This can be observed in 
Figure 4, which shows the situation for the gauge couplings 3 of the (MS)SM, g\ = 0.35, 
#2 = 0.62 and gj, = 1.23, for U(1)y, SU(2)l and SU(3) C interactions, respectively. 

Of course the variance of L or H is not exactly the same as the variance of the flat 
direction, but still it gives a reliable measure about the freedom of the flat directions 
to fluctuate. To completely characterize the statistical properties of the flat directions, 
we should rather obtain instead the 2-dimensional PDF V(L,H), which describes the 
simultaneous distribution of both H and L. In particular, V(L' , H')dLdH represents the 
probability of finding (L,H) within the field-space region [L',L' + dL\ x [H',H' + dH\. 
Therefore, similarly to the one-dimensional case, the V{L,H) function can be obtained 
from the numerics, simply by dividing the 2-dimensional unit-normalized histogrammatic 
distribution of the fields, at each time by the field bin steps AH and AL. 

3 To be more precise, one should let the gauge couplings run to the scale of inflation ~ -J HM P . However, 
for our purposes, this would not change - qualitatively speaking - anything; the U(l) and St/ (2) couplings 
would simply grow slightly, and thus the flat direction fluctuations in that case would result slightly more 
blocked, whereas although the SU(3) coupling would indeed decrease, it would never become small enough 
for the flat direction to fluctuate as if it were a free massless field. 
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Figure 5: Sequence of snapshots at different moments of the evolution of the 2-dimensional PDF 
V(L,H). From left to right we consider the case of having a U(1)y, SU{2)l and SU(i) c gauge 
coupling. From top to bottom, the rows correspond to N — 5,40, 150 and 1000 efolds (note the 
change of scale in the field space for different efolds) . Looking at the snapshots in the bottom row, 
one can appreciate how little the flat direction fields fluctuate at large couplings. 



In Figure 5 we show the evolution of the field fluctuations in the (L, H) plane for all the 
(MS)SM gauge couplings, by plotting V(L,H) at different moments. One can clearly see 
the probability leaking into the flat directions L = ±H. As expected from the behaviour 
of the variances shown in Figure 4, such leaking probability is of course higher the weaker 



- 13 - 



Prob \ N 


AT = 40 N = 150 iV = 1000 


Prob(|# ±L\ > 2;g < 1) 
Prob(|iT±L| >2;g u{1) ) 

ID—— V./'l it i r 1 v n, „ \ 

Probdif ± L| > 2;3 St/(2 )J 
Prob(|if±L| >2;g su{3) ) 


14 % 59 % 92% 
0.30 % 8.2 % 29 % 

0.2 r 70 0.0 70 21 70 

0.16 % 3.8 % 9 % 


Prob(|#±L| > 5;g <C 1) 
Prob(|tf ±L| >5; 5t/(1) ) 
Prob(|tf ±L| > 5; 5sc7(2) ) 
Prob(|tf ±L| > 5; 3sc/(3) ) 


0.04 % 12 % 73 % 
0.009 % 0.56 % 12 % 
0.0010 % 0.35 % 7.4 % 
0.0007 % 0.11 % 1.1 % 



Table 1: Here Prob(|iJ ± L\ > A;g c ) represents the probability of finding any flat direction with 
an amplitude with absolute value greater than AH a , when the coupling is g c . The case j « 1 
represent the free field case. 

the coupling. This is visualized in Figure 5, where we have fixed the scale in the horizontal 
to be same in all the snapshots taken at a given time for each different couplings. One can 
very clearly appreciate that the 'branches' developing in V(L, H) along the flat directions, 
reach greater distances in the (L, H) plane the weaker the gauge coupling is. For the U(1)y 
case, the weakest of the (MS)SM couplings, the flat direction branches reach values as high 
as 12H . But even in this case, in which the flat directions are the least constrained to 
fluctuate (compared to the SU(2)l and SU(3) C cases), such high amplitude regions are 
still highly unlikely. The 2-dim PDF obtained for the different coupling regimens are very 
suppressed for those high field values, if we compare them with the PDF for the case of a 
massless flat direction. 

To assess the degree of blockage of the flat directions in the different coupling regimes, 
we show in Table 1 the probability of the fields to be above certain fixed thresholds. We 
then compare with the probability that one would expect for a massless non-interacting 
flat direction field. The numbers speak for themselves. 



5. Conclusions 

During inflation, all the scalar fields fluctuate. The non-flat directions are coupled with 
each other either by virtue of Yukawa interactions or, as is often the case, via the D-term. 
Their variances will spread out to values of the order of a fraction of H Q in a matter of 
relatively few efolds. These variances act as effective mass terms for the flat direction 
(f>, whose fluctuations therefore become blocked and eventually saturate as equilibrium is 
reached with ((p 2 ) = aH 2 , a ~ 0(1) — 0(10), unless the coupling is very weak. The stronger 
the coupling strength, the more effective is the blocking, as is evident in Figures 3 and 4. 

We have verified these expectations quantitatively in Sect. 3, where we set up the 
the coupled system of Langevin equations, which account for the field fluctuations during 
inflation, and studied the fluctuations numerically. For simplicity, we focused on the F- 
terms and D-terms separately. Our main results are encoded in Figures 3 and 4, which 
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elucidate the spreading of the field variances as the function of the number of efolds for 
different coupling strengths. 

We have assumed that initially, at the beginning of inflation, all the fields are at the 
origin of the potential. Any random displacements around the origin are not likely to 
change qualitatively our conclusions, except of course for those initial conditions where the 
flat direction field has from the outset an amplitude much larger than H Q . It should also be 
emphasized that we assume inflation in the background with all the susy scalars acting as 
spectators. This is not the case in MSSM inflation [14], which makes use of a saddle point 
along the flat direction. Thus, the present considerations do not apply to MSSM inflation, 
where the initial amplitude is put in by hand. 

However, many other cosmological deliberations that tacitly assume a large amplitude 
along the flat direction generated during inflation should be reconsidered. There are two 
issues one should bear in mind when discussing flat directions and inflation: the blocking 
of fluctuations due to non-flat directions, discussed in the present paper, and the number of 
efolds actually needed to reach equilibrium, which concerns the inflationary sector rather 
than the more precisely defined susy scalars. 

Finally, although we have considered F-flatness and D-flatness separately, it might be 
possible to solve the Langevin equations for the full MSSM scalar field contents and follow 
the fluctuations of all the fields, including the simultaneous fluctuations along many flat 
directions. This would be an interesting exercise. 
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Appendix. Some aspects of the derivation of the Langevin equations 

Let us consider some physical scale L and IR/UV decompose any field (j) with respect to 
the physical volume Vl ~ 1/L 3 . Choosing our frame of reference as centered in such a 
region, the IR part <fi m is coordinate-independent within Vl, and it can be identified in 
fact with the volume average of (j> over Vl, i.e. 



This holds as long as the window function W used in the IR/UV decomposition has compact 
support in real space. In de Sitter space, the natural choice is L = H' 1 . Considering a 
region of physical volume Vh ~ each field can then be decomposed as (fii(x,t) = 

$>i(t)+ipi(x., t) , with <E>j and fa obtained according to eqs. (3.1)-(3.2), as long as |x| <C H^ 1 . 




(5.1) 
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The IR part <I>j can be indeed interpreted as the contribution to fa from all field wavelengths 
larger than H~ l at the given region. 

Let us now derive the Langevin equations in some detail. First of all, note that the 
IR/UV parts of the derivative of the field are not the same as the derivative of the IR/UV 
parts of the field, i.e. (IR{ cp})' (x, t) / IR{0} and (UV{0}) (x, t) / UV{</>}. This curious 
property, which emerges simply from the fact that in an expanding universe the window 
function W(k, t) for a fixed physical volume depends on time, has nonetheless profound 
consequences for the dynamics of fields living in De Sitter space. This property is indeed the 
key factor behind the derivation of the Langevin eqs. describing the stochastic behaviour 
of the IR dof. 

In the hamiltonian picture the eom are 



7Ti , TTi + 3H TTi = -^V z <pi - diV , (5.2) 



with diV = dV/d(j)i. We should thus IR/UV decompose independently 4>i and m as 
4>i(x,r]) = <&i(rj) + </?j(x, 77) and 7Tj(x, 77) = IT (77) + <57Tj(x, 77), and introduce such decomposi- 
tion into eqs. (5.2). Assuming that V({(frj}) is infinitely differentiable with respect to any 
field 4>i, then 



diVMj}) = DiV({^}) + £(AjV>i + • • • , (5.3) 



with Di = d/d$i, Dij = d 2 /d<S>id<£>j, and V = V({$j}). We should then plug this 
expression for diV into the eqs. (5.2), and separate the evolution of the IR dof from the 
UV dof . One finds this way 

/<i 3 k 
^pe-* kx <MM^(M) (5.4) 

1 - \ - f <i 3 k 

Ui + 5m + 3H {Ui + 6m) - ^VVi + DtV + X)(AjV>j + - = J e~ llsx m{^,t) W(k,t) , 

(5.5) 

where we have neglected V 2< l?j, since this term is sub-dominant in determining the evolution 
of the average value of the field <l?j ~ (4>)v L inside a volume Vl ~ Ignoring short 

wavelength interactions, the UV dofs satisfy 

(pi - 6m = (5.6) 
5m + 3H a 5m - \v 2 ifi + DiV + ^(DijVfaj + • • • = , (5.7) 



4 If we consider, however, the derivative of the total (f>, and not the derivative of each of its IR and UV 
parts separately, then it is true that 0(x, t) = (<^ IR )'(x, t) + (<^> uv )'(x, i) = IR (x, t) + <f) uv (x, i). 
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so, subtracting these eqs. to eqs. (5.4)-(5.5), we are then left with the equations for the IR 
dof as 

r j 3 w 

$i = Ik + J ^3 e" 4kx &(k, t) W(k, t) (5.8) 

H* = 3H U t + DiV + J e~ lkx 7Ti(k, t) W(k, t) (5.9) 

Considering the previous equations as operator equations, we arrive at the usual in- 
terpretation of the eq. of the IR modes as Langevin eqs (see main text). To completely 
characterize these equations, we need however to choose a window function. The simplest 
choice for it is a step function 

W(k, t) = 0(e aH Q - k) = 1 - 9{k - e aH ) , (5.10) 

where the number e < 1. Although this window function does not have a compact support 
in real space it suffices for our purposes as we will not be interested in the exact spatial 
correlations of In such a case W(k,t) = ea(t)H^5(eaH — k). Using this and the usual 
normalization of the creation/annihilation operators, 

<a t (q)4(q')> = (27r) 3 <5y$ 3 (q - q') , (5.11) 

(^(q^q')) = (4(q)a5(q')) = (aj(q)6j(q , )> = , (5.12) 

one arrives at 

sf (x,x',M') = ^a\t)\f4 q ,t)\l €aH J (aH\^-^\)5(t-t'), (5.13) 

sJjWiM') = ^aHt)\U( q ,t)\l taH j (aH\^-^\)5(t-t'), (5.14) 

where jo(x) is a spherical Bessel function of order 0. Had we chosen another function 
for W(k,t), like a gaussian or some other function smoothed around the scale eaH, then 
correspondingly we would have found, instead of a 'sharp function' 5(t — t'), rather a smooth 
function peaked at t = t' . Moreover, note a couple of things. First, we only need and 

sjj at x = x', since we are only interested in the random kicks that an IR mode $j receives 
in its domain of size ~ H~ s , and not in the correlation between kicks at different (causally 
disconnected) domains. Therefore, the spherical Bessel functions will be dropped, since 
jo(0) = 1. Secondly, these correlators are indeed divergent at equal times t = t 1 , as a 
reflection of the fact that we have chosen a 'sharp' step-function for W(k,t). However, the 
presence of 5(t — t') will not be a problem for our purposes, since as explained in Section 
III, the correlators we really need to compute are 

/t+dt pt'+dt tt4 pt+dt 

J s\f(t,t')dtdt' = e%^ 2 J dta 3 (t)\U( q ,t)\l faHo , 

(5.15) 

rt+dt pt'+dt TTi pt+dt 

S%\t,dt) = J J s ( i ?(t,t>)didt> = e% j ^j t dta 3 (t)\U(qM q=eaHo , 

(5.16) 
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In a pure de Sitter background for Inflation, the solution to the mode equations (with 
boundary conditions matching Minkowski modes at k — > oo) are [13] 



d (x 3 / 2 H^{x) 



2 fc3/2 r ' Jny ' 1 ~ 2 kVi a (t) dx 



(5.17) 



with H^X(x) a second- kind Hankel function of order |, x = kr\ and Tj = J dt/a(t) the usual 
conformal time. Using these mode functions evaluated at k = eaH Q , one can therefore 
obtain the exact form of the correlators S;f and Sffl. In the limit in which dt <C 1/H , 
one finds 

s\f{t, dt) = (1 + e 3 )5 tJ ^dt , dt) = e% 3 j0^dt , (5.18) 

We have used these noise correlators evaluated at e = 0.1, and observed that the results 
do not change if one makes e even smaller. 
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